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Abstract 

We present a systematic way to compute the scaling exponents of the structure functions of the Kraichnan 
model of turbulent advection in a series of powers of ^, adimensional coupling constant measuring the degree of 
roughness of the advecting velocity field. We also investigate the relation between standard and renormalization 
group improved perturbation theory. The aim is to shed light on the relation between renormalization group methods 
and the statistical conservation laws of the Kraichnan model, also known as zero modes. 
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1 Introduction 



Fully developed turbulence seems to have properties that are familiar from an another branch of physics, the theory 
of critical phenomena. Thus certain observables are scale invariant, i.e. they exhibit power-law dependence on 
length scale. This power-law seems to be reasonably well captured by dimensional analysis, with however systematic 
discrepancies occurring that don't seem to have simple structure. This resembles the phenomena of second order 
phase transitions where dimensional arguments (mean field theory) do a reasonable job, but don't fully account for 
the true scaling exponents. Furthermore, in both cases the observed scaling exponents seem to exhibit universality, i.e. 
a relative independence on many details of the system: in the case of critical phenomena details of the microscopic 
Hamiltonian are unimportant, only symmetries matter, in the case of turbulence, microscopic details of the forcing 
mechanism that maintains the turbulent state seem irrelevant. 

Important differences occur too. On the practical level, in the theory of critical phenomena a rather explicit starting 
point for calculating correlation functions exists in the form of Gibbs measure given in terms of an explicit Hamilto- 
nian. This is directly accessible numerically and analytically one can gain qualitative understanding by perturbative 
study around the upper (and sometimes lower) critical dimension. In case of turbulence, while the dynamical equa- 
tions governing fluid motion have been known for long, the analog of Gibbs distribution is not. The nature of the 
stationary state describing temporal averages of measurements is a dynamical problem that is unsolved. Also, there 
doesn't seem to be a parameter in the problem whose special value would make the dimensional analysis exact (like 
the dimension of space in critical phenomena) and which could then provide a basis for a perturbative study of the 
problem of anomalous scaling. 

On the more fundamental level there also are differences. The modern renormalization group (RG) theory of 
critical phenomena is based on locality in position space. The effective theory of any given scale is given by a Gibbs 
state defined by a local Hamiltonian. The RG relates these different effective theories to each other. In turbulence 
the stationary state is characterized by fluxes ("cascades") of conserved quantities (energy and in 2d enstrophy). This 
cascade process is believed to be local in wave number space and it is not clear what the right RG description is. Both 
direct 1251 1261 1211 1501 and inverse (see for review and also |4l for criticism) RG's have been proposed in the past. 
The former is analogously to the theory of critical phenomena based on locality and coarse graining in physical space 
whereas the latter coarse grains in wave vector space preserving locality there. 

Theoretical progress in critical phenomena came in two ways: by exactly solvable models with nontrivial scaling 
(the 2d Ising model) and by perturbative RG analysis near a Gaussian theory (e-expansion). For Navier-Stokes turbu- 
lence these options are not available. The analog of the e-expansion is provided by considering random stirring where 
the power spectrum of the force concentrates at large wave numbers k, being proportional to A;^"'^^^'^. Unfortunately 
for small e this is very different from the turbulent situation where the force is concentrated in low wave numbers. 
Nevertheless interesting lessons can be learned also from such small e-expansions as witnessed by the renormalization 
group studies carried out by the St. Petersburg school (see for example |4l|M| and references therein). 

A different approach which can be pursued in the investigation of fully developed turbulence is to consider phe- 
nomenological models able to capture some of the properties of turbulent fluid motion. Among such models, the 
Kraichnan model of passive advection l35l l36ll has permitted in the last years to shed some light the mechanisms 
underlying the genesis of intermittency l30llT7lfT0lfT8l[mi28ll . An exhaustive review of the results derived from the 
Kraichnan model can be found in |24]. 

The simplifying assumption which defines the Kraichnan model is the replacement of the Navier-Stokes velocity 
field advecting the scalar observable with a random field, Gaussian and (5-correlated in time. This latter assumption is 
crucial as it ensures that the correlation functions of the scalar field satisfy closed Hopf equations that allow to relate 
the n-point functions to n — 2-point functions via Green functions of differential operators built out of the spatial 
part of the velocity covariance. The latter involves a parameter ^ describing the smoothness of the velocity field: its 
realizations are Holder continuous with exponent less than ^/2. 

The properties of the theory versus ^ are particularly interesting. At ^ equal zero, the effect of the advection 
is to "renormalize" the microscopic molecular diffusivity to a macroscopic eddy-diffusivity. The resulting theory is 
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Gaussian and provides a natural starting point for a perturbative investigation of the system in powers of ^. Although 
it can be argued that the value of ^ ideally corresponding to a turbulent flow is equal to 4/3, the scalar field tends 
already for small values of ^ to a steady state where an inertial range sets in. Thus, the turbulent regime it is accessible 
in perturbation theory. This is at variance with what happens for the e-expansion of the Navier-Stokes equation where 
the perturbative expansion has its starting point in a model with vanishing inertial range. 

The scalar correlation functions were seen to have the zero molecular diffusivity limit order by order of the 
^-expansion. This result was subsequently proved rigorously for all ^ in ITTIl . However, the main result of the ^- 
perturbation theory has been the derivation of corrections to the naive scaling prediction of scaling exponents of the 
structure functions of the scalar 1 30. . 1 0. .32. . 1 1 . .28. .29 1 a result that also was obtained in a perturbative expansion of 
the structure functions in inverse powers of the spatial dimension IITtI ITS 44 |. 

Both perturbative approaches were based on the study of the Hopf equations satisfied by the equal time correlations 
of the scalar. In the inertial range where forcing and dissipative effects are negligible, these equations reduce in the 
stationary state to the annihilation of the correlation function of order n by a linear operator For each n, the 
operators A1* admit scaling zero modes that can be computed perturbatively in ^ (or in 1/d ) and can be shown to 
determine the leading scaling behavior of the structure functions. 

The zero modes are statistically conserved quantities of the scalar field. The presence of such conservation laws 
provides a mechanism underlying the phenomenon of intermittency. Indeed, they were observed numerically both in 
passive scalar advection by a two-dimensional Navier-Stokes velocity field [16J and in shell models [58 6 9| (see 
also EUCl). 

The concept of zero mode has therefore proved fruitful both to shed light on properties of more realistic models of 
turbulence and in the analysis of the Kraichnan model. However, an issue which was left open is how to compare the 
small ^ expansion of |30, 10 1 with the field-theoretic methods based on perturbative expansions of the Martin-Siggia- 
Rose functional that describes the time space time correlation functions of the theory. The interest of such question is 
threefold. 

First, the Martin-Siggia-Rose formalism allows to study all the correlation functions described by the theory 
without restriction to equal times. Hence it permits to inquire what the standard perturbation theory in ^ can say about 
more general observables. 

The second reason of interest is that the Martin-Siggia-Rose provides a natural framework to define an algorithm 
to compute, at least in principle, all higher orders of the perturbative expansion. In the zero mode approach, scaling 
exponents have been computed using scaling Ansatze, which have provided the first order in ^ term (or in some 
special cases exact results, see for example l52l ). On the contrary, using the Martin-Siggia-Rose formalism, ultra- 
violet renormalization and short distance expansion, expressions up to 0(^^) of the scaling exponents were derived 
in |i5||2|(3]J. 

The fact that ultra-violet renormalization proves useful in the context of the Kraichnan model might appear at first 
glance surprising. In general, renormalization is used to make sense of perturbative expansions affected by ultra-violet 
divergences. On the other hand the Kraichnan model, at least as it should be defined whenever the small ^ expan- 
sion is used, exhibits a well defined small scale behavior. It is therefore interesting to establish a precise connection 
between the zero-mode and the ultra-violet renormalization methods. Zero modes provide the leading inertial range 
asymptotics of equal time correlations. It was pointed out already in flTIl that taking the Mellin transform yields a 
well defined way to identify the zero mode contribution given the full, all scales, expression of correlation functions. 
In the present paper this idea is developed in order to show how it can be implemented in principle to all orders in 
perturbation theory. As an example, explicit expressions of isotropic and anisotropic scaling exponents of the structure 
functions are obtained to second order in ^. The result does not just recover the result found with the renormalization 
group in |5l but establishes the relation between the two methods and the reason why the predictions must coincide to 
all orders in perturbation theory. 

The third reason of interest is the relationship between the different renormalization groups, direct and inverse, 
that have been proposed for turbulence. In particular the concept of inverse or infra-red renormalization has been 
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proposed as an alternative tool to solve scaling in the Kraichnan model l3ni32ll . The special features of the Kraichnan 
model are probably not suited to settle the issue which of the approaches is more natural. However it provides a simple 
case study of what infra-red renormalization is about. 

The scope of the present paper is to address in a comprehensive way the these three issues. 

The paper is organized as follows. In the first section|2lthe passive scalar equation is introduced and the definition 
of the Kraichnan ensemble is given. In section |3l the asymptotic expression of the velocity field in the physically 
relevant ranges are derived. The analysis of the correlation of the velocity field shows that the roughness parameter ^ 
lend itself as a parameter for a perturbative construction of the solution of the Kraichnan model. 

In section|4]the Hopf equations governing the dynamics of equal time coiTclation functions are recalled. Section|5] 
delves into the relation between the Mellin transform of the solution and statistical conservation laws. Statistical 
conservation laws are specified by homogeneous solutions of the Hopf equations admissible in the inertial range but 
not matching the boundary conditions at large spatial scales. Such solutions, referred as zero modes can be interpreted 
as the residues of the first poles of the Mellin transform of the full solution of the Hopf equation. 

In section|6lwe recall the analysis of the zero modes in (TTi which allows to predict the inertial range asymptotics 
of the structure functions of the scalar field. The prediction is encoded in a scaling Ansatz which can be tested in a 
perturbative expansion in powers of the Holder exponent of the velocity field. 

The perturbative expansion is couched in the field theoretic formalism through the introduction of a Martin-Siggia- 
Rose generating function (section^. The limit of zero Holder exponent is shown to provide a Gaussian limit around 
which it is possible to develop a perturbative expansion according to standard Feynman rules. 

In section [8l the explicit expression of the scaling exponent is given up to second order in ^. The calculation 
is based on the use of Mellin transform technique described in section |6l Technical details of the evaluation of the 
integrals are deferred to appendices ICl iDland |53|. 

A consequence of the scaling Ansatz of section|6lis that the scaling exponent of the structure functions also govern 
the blow up rate of scalar gradients versus the dissipative scale. This is recalled in section |9l where the exponents are 
again evaluated to the second order in ^, the evaluation being considerably easier than via the structure functions. 

In section no|the Wilson's formulation of the renormalization group is recalled first in the traditional direct form 
and then in the inverse setup. In section^Jthe infrared scaling fields exhibiting the anomalous scaling of the structure 
functions are derived and in section EJthe ultraviolet scaling fields responsible to the scaling of the scalar gradients. 

The last section is devoted to conclusions. In an appendix some details of computations are collected. 

2 Passive advection by a turbulent fluid and the Kraichnan model 

The passive scalar equation describes a scalar quantity 0(x, t) which is advected by a fluid moving with velocity 
v(x, t) and diffuses with molecular diffusivity k: 



The molecular diffusivity k describes microscopic dissipative effects. The role of the force term / is to provide a 
source for the scalar in order to sustain the system which otherwise would decay under the effect of dissipation. 

A turbulent steady state with large inertial range sets in if the forcing acts over spatial scales much larger than 
those where dissipation becomes relevant. The details of the forcing are supposed to be irrelevant for inertial range 
scaling. It is therefore convenient to chose the forcing as a Gaussian field of zero average and covariance 



(1) 




(2) 



The spatial part of the covariance F is a smooth function decaying rapidly at infinity and satisfying 



F (0) = > 



(3) 
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Hence the constant L p specifies tlie cliaracteristic scale (tiie integral scale or tlie correlation length) of the forcing. 
The delta correlation in time in ^ means that Q is a stochastic (partial) differential equation. 

In realistic models of turbulence the velocity field is specified by a solution of the incompressible Navier-Stokes 
equation. One would like to understand the typical behavior of the scalar given a typical realization of v from an 
ensemble of such solutions. In the Kraichnan model |36| this ensemble is replaced by an ensemble of Gaussian 
random velocity fields which is chosen so as to mimic some properties that are thought to be crucial of real turbulent 
velocity ensembles. 

One considers random velocity fields with Gaussian statistics having zero average and covariance 

^ v"(xi,ti)v'3(x2,t2) ^=5(t2-ti)I?"^(xi-X2;m,M). (4) 

Real turbulent velocities have two characteristic length scales, the short, dissipative scale M^^ and the long integral 
scale m^^. These are modeled in the Kraichnan ensemble by an ultraviolet cutoff M and infrared cutoff m in wave 
numbers entering the spatial part of the covariance in @. Moreover, realistic turbulent velocities are approximately 
self similar for scales between the dissipative and the integral scales. 

The simplest choice for the spatial part of the covariance having such properties is 

^-^u''/^(cDxi^^M]{q) (5) 

where 

n"^(q) :=r^-q"q'^ (6) 

with q denoting the unit vector q_/q, q = |q| and X[m,M] il) is the characteristic function of the interval [m, M]. 

We consider in this paper only incompressible velocities 5 • v = which is guaranteed by the tensor 11°^. It is 
convenient to introduce the covariance of velocity differences: 

-< K(0,t')-v°(x,t')][v^(0,t)-v'^(x,t)] ^=2 5(t-t')rf"^(x;m,M) (7) 

where we have introduced the spatial velocity difference covariance 

d"^(x; m, M) = D^f (0; m, M) - D^" (x; m, M) = Dq^ J ^ n-^(q)X[™,M] (q) (8) 

As seen in detail later, d"'^ is approximatively scale invariant in the inertial range: 

d°^(Ax; m, M) ~ A« d°^(x; m, M) (9) 



for 



The constant 



where 



M ^ < |x| ,A|x|, A|x| ,|x| < m ^ (10) 
D"^ {0;m,M) = D{m-^ - M-^)d''^ (11) 



^-^°^(^' ^'-rm ^^^^ 

describes the mean square velocity field which blows up as the integral scale m^^ tends to infinity. 
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Finally, the delta correlation in time of the velocity fields guarantees the statistical invariance of the velocity 
differences under Galilean transformations v'(x, t) = v(x + ut, t) — u, an important property of the Navier-Stokes 
equation. More important, it leads to a relatively explicit solution of the scalar statistics. 

Equation © together with Q defines an infinite dimensional stochastic differential equation with multiplicative 
noise. In order for this object to be well defined according to the general rules of stochastic calculus l47ii22J . it is 
necessary to specify how the product between the velocity and the scalar field in Q is defined as the continuum limit 
of a finite difference stochastic equation. The Kraichnan model is supposed to be the limit of a physical system with 
finite time correlations. Hence it is natural to regard ©, @ as infinite dimensional stochastic differential equation in 
the Stratonovich sense. 



3 Asymptotics expressions of the velocity ensemble 

Far from being the only possible, the choice ^ suits the derivation of explicit asymptotic expressions of the velocity 
field. Of physical relevance is the behavior of the velocity field in the dissipative and inertial range. 

3.1 Dissipative range asymptotics 

The dissipative range is defined by the inequalities 

mx <C Mx < 1 (13) 

Under this assumption the Fourier exponential in ^ can be expanded in Taylor series. Up to leading order, the 
expansion yields 

having introduced the rank-two real-space tensor 

r"'3(x,z) =5"^-— ^ x"x^ (15) 

a — 1 + z 

Whenever 

m 

- « 1 (16) 

the eddy dijfusivity 

x^:=Dm-^ (17) 

dominates the velocity field in the dissipative range. The leading correction to the constant mode of the velocity field 
is smooth and vanishing as M^^ tends to zero. 

3.2 Inertial range asymptotics 

The integral ^ is convergent both if the ultra-violet cut-off M is set to infinity and the infra-red cutoff m to zero. 
Thus, the statistics of the velocity differences exists in such limit. It can be determined by considering the Mellin 
transform of the spatial correlation of the velocity field 

r°° rim 1 C rl'^n p*"' q-x 



= 0„{ I -- J ^-^n-«i) (18) 

q>Tn 
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The integral is convergent for 5?z < and can be performed explicitly. The result is (see appendix |B] for details): 

Z—i; Z 

D''P{x;m,z) = D^ciz) '^ _^ T°^(x,z) (19) 

z [z — t,) 



with D and T" ^ respectively defined by (fT2l and (fTSt whilst 

( d-i + z)r(^)r(i 
(d-i) 2-r(^^) 



is a function with simple poles for z a positive even integer . 

The small scale asymptotics is derived by evaluating the inverse Mellin transform involving an integral over z 
along = const < by pushing the contour to the right and picking residues from the poles. This gives 

D"^(x;m) =Z)m-«5"^-(i*"^(x)+o(m2-«a;2) (21) 

Thus, the residue of the pole at zero corresponds to the eddy diffusivity. The pole at ^ specifies instead the inertial 
range asymptotic of the structure tensor of the velocity field: 

d*"/3(x) := D c(e) x« T"'3(x, (22) 

Note that the blowup of c(^) at ^ = 2 is related to the fact that if ^ is equal two the pole at z equal ^ in ( fT9l turns 
from simple to double. This indicates the existence of logarithmic corrections to the analytic behavior, proportional 
to x^, of the velocity field structure function. Logarithmic corrections are suppressed for example by redefining 

D„ ^ (23) 

r(i-i 

The rescaling does not affect universal quantities in the small ^ limit and will be neglected in the present paper. 

Let us finally discuss the behavior of the velocity covariance around ^ = 0. At fixed ultra-violet cut-off it vanishes 
linearly with ^. On the other hand, the removal of the ultra-violet cut-off reduces the dissipative range to the the single 
point X equal zero. There the velocity field coincides with the large scale constant mode. In the inertial range, by dlTTl 
the velocity field is vanishing with ^ for any nonzero m and x also after sending M to infinity. We will see that the 
expansion around ^ equal zero provides a viable analytic tool for the investigation of universal properties of advection. 



4 Hopf 's equations and statistical conservation laws 



Let 0{-K,t) be the solution of the stochastic differential equation Q with suitable initial condition. The equal time 
correlation functions: 



2n 



C2n :=-< n^(^^'*) ^ 



(24) 



i=l 



satisfy in the Kraichnan model a solvable hierarchy of Hopf equations which are simplest to derive in the Ito repre- 
sentation BTIE^ of the equation. Letting v • d^O be defined with the Ito convention the equation Q becomes 



1 



By (fm . in the Ito representation the molecular viscosity is renormalized by the velocity field to 

K= K + D {m^^ - M~^) 



(25) 



(26) 
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In the inertial range where dissipative effects are neghgible the diffusion felt by the scalar field reduces to the eddy 
diffusivity x^, introduced in equation dTTl . 

A direct application of the Ito formula then yields 

dtC2n -^Yl C^Sn - Yl ^"^(^^ - ^)d^fd^^C2n = T^n (27) 

^2n ■■= C(2„-2) ^ ^ = E C(2n-2)(X1, ■.■,^2n,t) F (^^^^) (28a) 



with 



T^:=f( ^^V^ ) (28b) 



Lf 

Odd order correlation functions will be ignored since they vanish in the steady state due to parity invariance 
For a translation invariant initial condition for 6 (say 0) C2n is translation invariant i.e. 



^5x^C2„ = (29) 



the Hopf equations reduce to the final form 



dtC2n -^y\ dl<Z2n + V d"^(x, - X,; m, M)d^<^d^,C2n = ^2n (30) 

t Kj 

Note that owing to translational invariance the 2n-point correlation in d dimensions is a function of dn = (2n — l)d 
variables. 

Since d30b depends on the velocity statistics only through the structure function it has coefficients that have limits 
as the ultra-violet and infra-red cut-offs of the velocity field are removed. In that limit we may define the elliptic 
negative definite operators 

M'2n :=5^d*"^(x.-x,)ax^9^, (31) 

i<j 

Since the coefficients d* ^ are homogeneous functions of degree ,f these differential operators are self similar, namely 
homogeneous of degree ^ — 2. The stationary state correlation functions of the scalar satisfy the equations 

2n 

- ■^2nC2n = ^ Yl C^2n + ^2 n (32) 
^ i=l 

The terms collected on the right hand side of d32t are related to the non-universal dependencies of the dynamics. 
Self-similarity breaking from the large scale can occur through the integral scale of the forcing Lp. Comparing the k 
and the Al2„ terms hints at the existence of a scalar field dissipative scale 

The existence and uniqueness of the solutions of (I32t has been proved at zero molecular dissipation for all values 
of the Holder exponent ^ < 2 of the velocity field L37.I . They are given as 



C2„(X;Li.) = - J dYM^2{X,Y)T2n-2{Y;LF) 



(34) 



7 



for X,Y G The kernel of the operator Al2n ^ proved to be locally integrable and the integrals 

converge absolutely as long as the forcing scale is finite. By the homogeneity of A^2n~^ S^*^ immediately that 

C2„(X,Li.) = L'f-^^C2n {^^, l) (35) 
i.e. the canonical dimension of C2n is ra(2 — ^). 

5 Zero modes and short distance asymptotics 

From d35l we see that the large Lp behavior of the correlation functions C2n is dominated by the large scale velocity: 
they blow i 
transforms 



they blow up like L^J""^ ^''C2n(0, 1). Sub-leading terms in the inertial range can be extracted by considering the Mellin 



C2n{^,Lp;z):= = / — --^^"'-^ (36) 

w Lip /n w 



These integrals converge for the real part of z small enough and are expected to extend to meromorphic functions of 
z at least for generic ^. 

The Mellin transform of the Hopf equation ( 1^ is 

- M*2nX'+^-^C2n{±, Lp; z + 2 - ^ = |5iX^+2j^„(X, Lf- z + 2) + X'p2n{X Lf; z) (37) 

It was observed in lITTl that poles of C2n can occur either for values of z for which J^2n has a pole or for z such that 
the operator M.2n h^s a zero mode 

A^2n3 = (38) 

which is a homogeneous function of degree z. The poles of ^2n are in view of (I28at determined by solving the Hopf 
equations lower in the hierarchy. One then ends up with an asymptotic short distance expansion for C2 n in terms of 
homogeneous functions of the coordinates 

C2„(X, Lf) = X^-' L''/-^^-'^"" Aj{X) (39) 
j 

Since the forcing covariance is smooth and has a Taylor expansion F{x) = Yl /n|3;|" we conclude that the scaling 
exponents (nj that may enter in ( l39l are either the homogeneity degrees of the zero modes of A^2n '^^ they are 
determined in terms of the previous ones Cn-i,j- 

The general situation is illustrated by the case of the two-point function C2 which may be computed explicitly by 
quadrature if the forcing is isotropic LIU : 



~ /o"^ dx2F{x2/LF)xt ' 

xf-^ + + k] 



C2{x;Lf) = + 6 / dxi , / T Z \, , (40) 



with (i*"^ specified by In such a case, at zero molecular viscosity, equation ( BTl reduces to 

- d*'"^{jc)da,df3x'+^'^C2{l, Lf;z + 2-0 = x' Lf; z) (41) 

with solution 



C2[x,Lf;z)--Lp —-^^-^-^—^^^-——^ 



(42) 



Zero modes correspond to the poles in z equal zero and 2 — d — respectively associated to short and large distance 
asymptotics [18 24 1. The only zero mode contributing in short distances is the constant, the poles of F are at 
z = 2 — + n for nonnegative integer n. 
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6 Anomalous scaling of the structure functions 



Correlation functions that probe the sub-leading terms in the short distance expansion (I39t are provided by the struc- 
ture functions of the scalar namely 

S2n{^,LF,i) [9{x)-9{0)f^ y (43) 

where we denoted explicitly the dependence on the forcing scale and k via d33l . Let us consider for simplicity 
isotropic forcing so that S2n is only a function of x = |x|. We may scale it out as 

52„(x,L^,^) = x(2-0n52„ (l,^,-) (44) 

\ X X J 

Suppose S2n{^:L,£) had a limit as £ tends to zero and L tends to infinity. Then we would conclude 

S2nix,LF,i) = (^l + ^ (^-^, , (45) 



i.e. the structure function scaling exponent would be n{2 — ^). This is the prediction of the Obukhov-Corssin P461 I20I 
theory, analogous to the Kolmogorov 41 theory [33 . 3H W7\ UJi for Navier-Stokes turbulence that predicts there a 
scaling exponent n/3 for the velocity n-point structure function. 

To discuss the validity of the Obukhov-Corssin theory, we note first that the limit £ — > exists for the correlation 
functions and thus for the structure functions by the results of |37|. The large Lp behavior depends on the nature of 
the terms entering the expansion ( l39l . i.e. the zero modes of the inertial operators. Structure functions are obtained 
from correlation functions by applying the finite increment operator Tx 

52n(x,t) = 2xC2„(xi, . . .,X2n,t) (46) 

The operator 2x = Yii '-x^ ^i^h ^ counting the number of fields in the correlation function generates finite field 
increments according to the rule ix/(y) = /(x) — /(O). From (l46l it follows immediately that the only zero modes 
of C2n contributing to S2n can be the irreducible ones, i.e. those depending on all the independent variables xi — Xj, 
i = 2 . . .n. 

The existence and the properties of zero modes were thoroughly investigated in l30lfT7lfT0lfT8lfmi24ll . It was 

shown in |30, 10 1 for small ^ and in 1 17 1 for large d that for each n > 1 there is a unique irreducible zero mode 32n- 
The irreducible zero mode mode has scaling dimension C,2n = (2 — ~ P^n with p2n > 0. In 1 1 L] further arguments 
were presented for the conclusion that this zero mode enters the expansion ( l39l and is the dominant term. This means 
that we may write 

where the function S2n has a nonzero limit as £ ^ and L^? — > oo: 

lim lims2„ ( := S2„(0 > 0. (48) 

Lp-^OOi^0 \Lp X J 

The ratio between the constants F^, and D, defined respectively in ^ and (fTTT l. is introduced for later convenience. 

It was also argued in OUllTni that the sub-leading terms for the asymptotics of 52 n have exponents well separated 
from the leading one for small ^, namely 



S2n{^,0,^]--=S2niO+oll^Y ] (49) 
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We may call this the scaling Ansatz for the Kraichnan model. 

It is useful to express this result in terms of the Mellin transform. Let us shift for convenience z by 2n which is 
the ^ equal zero theory scaling exponent: 



52„(x,L; z + 2n) := 



dw S2n{wx, L) 



Q W W 



z+2n 



(50) 



which converges for ^z small enough (actually < — 0(^)). No specific assumption is made on the origin of the 
integral scale L appearing in (1501 . The scaling Ansatz (1?^ then yields 



S2n{x,L; z) :-- 



(51a) 



with the function A having poles for values of z differing from o"2n by terms at least of the order 0{^^). Furthermore, 
A takes on the first pole of the Mellin transform the value 



(52) 



In the following sections structure functions will be computed using perturbation theory around ^ equal zero. It is 
useful to observe that the small ^ expansion of the Mellin transform generates a Laurent series in z. The residues of 
the expansion have simple relations with the expansion in powers of ^ of the scaling exponent: 



d_ 



ln{52„(x,L;z + 2n)L<} 



4n(0) 



(4n)'(0) , 4'n(0) 



+ 



+ 0(z°) 



(53) 



Analogous relations hold true to all orders in ^ and also for anisotropic forcing. Thus, the Mellin transform permits 
to extract, systematically at any order in ^, zero mode contributions. Furthermore, as it will be seen below, taking the 
Mellin transform of the Feynman diagrams generated by the Martin-Siggia-Rose formulation |41 1 of the Kraichnan 
model greatly simplifies the explicit evaluation of the corresponding integrals. The fact is well known in field theory 
as the Melhn transform has the effect to map a theory with massive propagators into a massless one L59..48J . 



7 Martin-Siggia-Rose formalism 

The perturbative analysis of the zero modes of the inertial operators -M2n "^be parameter ^ is based on the fact that 
the velocity structure function d22t becomes constant as ^ tends to zero and the operators become Laplacians. It can 
then be checked that the solutions of the Hopf equations are correlation functions of a Gaussian field 6. For small ^ 
the distribution of 9 should be given by a perturbation of this Gaussian. It is straightforward to derive a perturbation 
expansion from the Ito stochastic differential equation d25l . 

For this, let us observe that the solution of (125 1 can be written as 

0(x, t)= [ dS $K(x, t\y, to)9{y, to) + fds [ d% 9^(x, t|y, s)f{y, s) (54) 
J Jto J 

with 9^ solution of the stochastic differential equation 

[dt - ^dij $n = -V d^^ (55a) 
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with fH(x, t\y, t) = 6^'^^ (x — y) and the product is defined with Ito convention. This equation can be solved as a series 
in multiple stochastic integrals 



^it\s) = Y. •••/ Rit\tnMtn)dtn-dR{tn\tn-l)...yr{ti)dti-dR{h\s) (56) 

which converges in of the probability space of v (39||. Here R is the fundamental solution to the heat equation 



i?(x,t|x',t') = j |^e^^-(''-''')-t^^ (57) 

Inserting d56b to d54t and averaging over v and / leads to an expansion of the correlation functions in terms of 
integrals involving R and the co variances of v and /. 

A convenient way to generate this series is provided by the so called Martin-Siggia-Rose (MSR) formalism (see 
I15ll5ni59l and references therein). It provides a graphical representation of the perturbation theory analogous to the 
Feynman rules of quantum field theory. 

The MSR formalism is derived by introducing the generating function 

Z{3, j) =^ j m e^^ ^'H{dt9 + V • - ^dle - / - ij) y^j (58) 

where 

{j,9)= dt d'^xjO (59) 



and -< • >~yj denotes the average with respect to the velocity and forcing fields. By introducing of an auxiliary 
"ghost" field 9 the MSR functional becomes 

Z{J,3) =< j V[e~e] e^M,J}-A ^^^^ (60) 

with 

A = -t{e,{dt + ^-d-^^e)-t{ej) m 

Inverting the order of integration and averaging over the velocity and forcing fields leads to 

Z{j,j) =^ ^-M{e:9)-l{eM)+{jfi)+{B,3) (62) 

where 

Ai{e,e) = ]^ {eda9 , ^98^9) (63) 
and the average is with respect the Gaussian "measure" with covariance given by 



^ip-(x-x')~>tp^ \t-t'\ 

(27rp xp2 



-< ^(x,t)0(x',i') / 7:d^ —2 ^(P)- (64) 



and 

-< 9{x, t)9{^', t') >-G= iHoit - t')R{^, t|x', t') (65) 
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The free response function (l65l involves the "Ito-Heaviside function" 

which insures that equal time contractions in d63b will not occur as follows from the Ito convention used for the 
stochastic integrals in d56b . The generating function satisfies the normalization 

Z(0,0) = 1. (67) 

The perturbation expansion for Z{j, j) is obtained by expanding e"-^^ in powers of Ai and expressing the resulting 
expectations of ^ , ^ in graphical terms as explained in Section 8. Before going to that let us make two comments 
concerning the infrared cutoffs and the small parameter in the expansion. 

The free response function d65b depends on the eddy diffusivity (fTTl which diverges as the infrared cutoff m of the 
velocity field tends to zero. Thus this limit cannot be taken directly in d62b . In section |4]the properties of equal time 
scalar correlations were discussed in the case when m, = 0. This limit can be taken in the MSR formalism only for 
equal time correlation and structure functions. In particular it can be shown |45| that order by order in perturbation 
theory, if the integral scale of the forcing Li;' is kept fixed, the limit for m tending to zero exists and leads to the 
scaling predictions of \ 

On the other hand, in the MSR formalism it is natural to study the opposite limit when the integral scale of the 
velocity field is smaller than that of the forcing. As argued in section |6l anomalous scaling of the structure 
functions is a consequence of the existence of homogeneous zero modes of equation (l38l where the infrared cutoffs 
don't occur. Thus, it is plausible that the scaling exponents don't depend on the order in which the infrared cutoffs are 
removed. It would be interesting to spell this out more explicitly. In what follows we will keep m fixed and study the 
mp ^ ^ limit. 

Finally let us comment on the small parameter of the expansion. In section|3]it was shown that the spatial part of 
the velocity covariance vanishes almost everywhere (i.e. at x 7^ 0) as ^ ^ whilst being bounded from above by the 
eddy diffusivity x^. Thus apart from the constant mode the d63b can be viewed as a small perturbation in the limit ^ 
tending to zero. 

8 Small ^ expansion 

The general features of the small ^ expansion can be summarized as follows. The coefficients of the expansion 
are determined by integrals symbolically represented by Feynman diagrams. The basic ingredients of the Feynman 
diagrams are the free scalar correlation d64l and response functions (I65t together with the velocity correlation Q. 
They have the graphical representation 

< e{x, t)e(x', t') = • • (68a) 



-< 9ix, t)9i^', t') = • (68b) 

-< t;°(x, t)w'^(x', t')^ = . . (68c) 

where end-line dots represent external points. In order to exhibit the expansion parameter, Feynman diagrams are 
constructed by connecting free response and correlations lines through the 0{^^) part of the interaction vertex 



Ai = -Ax= \ (69) 



wv- 
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The bars transversal to scalar lines in d69l represent spatial derivatives. Finally, all contributions are reordered in 
powers of the Holder exponent of the velocity field by expanding the residual ^ dependence of the interaction. 

It is worth stressing that this perturbation series is ultraviolet finite i.e. the resulting integrals have a well defined 
limit as the ultra-violet cut-off M tends to infinity. Singular behavior is instead exhibited in the limit of infinite integral 
scale, m tending to zero. 

Structure functions are obtained by applying the finite increment operator d46b to the perturbative expressions 
of equal time scalar correlation functions. This operation removes order by order in perturbation theory all terms 
proportional to powers of the integral scale otherwise present for dimensional reasons in the correlation functions. The 
resulting perturbative expansion of structure functions contains only expressions exhibiting logarithmic divergences 
in the infra-red. 

In section |6l it was argued that S2n is in the inertial range a homogeneous function of the spatial separation of 
degree C,2n determined by the unique irreducible zero-mode of the Hopf equation of order 2 n. According to d53t C,2n 
can be straightforwardly evaluated order by order in ^ by taking the Mellin transform of the perturbative expression of 
S2n- Working under these assumptions, in the following two subsections the calculation of C,2n is outlined for the first 
two orders in ^. For simplicity we set M = oo in the calculation. More details are deferred to appendices ICl and IDl 



8.1 Zeroth order approximation 

The zeroth order of the perturbative expansion corresponds to neglecting the interaction term d63l in (I62t . In order to 
simplify the notation we keep the explicit non perturbative ^ dependence (fTTt in the eddy diffusivity x^. This trivial 
^ dependence can be always expanded a posteriori to check of the decoupling of the scaling exponents into the part 
predicted by canonical dimensional analysis and a part associated to self-similarity breaking in the inertial range. 

The second order structure function of the Gaussian theory is obtained from d64ll and In the limit Lp oo 
we get 



>f(x) 



lim lim 2 



(fp 1 



Lj^FiLpp) = ^ y (^(P • x)^ Fip) (70) 



(27r)'^ L>p2 

In particular, if the forcing correlation is isotropic we get 

5f(x) = 



dD 



Structure functions of higher order are given in the Gaussian limit as 



(2n)!, 
2"n! 



(71) 



(72) 



If the forcing is anisotropic, scaling properties of the structure functions are identified by expanding them in a func- 
tional basis invariant under rotations. In d-dimensions, this scope is achieved by resorting to an expansion in hyper- 
spherical harmonics I.54J . Under rather general conditions, the generic Gaussian structure function takes the form 



>fi(x) 



x'F^ 
dD 



Y,KjYjoi^) + 0{Lf) 



(73) 



with Kj some non-universal constants depending on the forcing. 



8.2 First order approximation 

To first order in ^ the structure fimctions require the evaluation of the diagrams 

V{i;2) = 2:x._/r\^ (74) 
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and 



V, 



(i;4) 



1. 



(75) 



In dT^ . ( 1751 the increment operator is given by eq. d?^ . As L^? ^ cxd the momentum transfer along free correlation 
lines is suppressed. In such a hmit the two diagrams can be recast in the form 



V, 



(0) 

(i;2) 



v: 



(O)a ^2^(0) 



(1;4) a 



(76a) 



,(0)a/3 



(0)^ 



:(0)> 



using the general notation 



(n). 



v: 



(76b) 



(77) 



to count the number of derivatives with respect to ^ of a diagram and 



v: 



a/3 

(i;4) 



2Dom^ f d<^p 1 - e*P-^n°^(p^ 



• — ^AAH— 



{2'kY 



X]m,oo] (P) 



(78) 



In dTSl the convention is adopted to denote with truncated scalar correlation lines in the diagrammatic representation 
the absence of momentum transfer along those lines. The factor two in dTSl stems from the action of the finite 
increment operator Ix on the two response functions in the integrand. From dTSb it is readily verified the expected the 
ultra-violet convergence of vf^.^ ^ as well as its logarithmic divergence in the infra-red. 

The sum of (I76at and (I76bt weighted by combinatorial factors yields the inertial range asymptotics of structure 
functions 



52n(x,m 

4 



(2n)! 



5f(x) 



in{n-l) reW. Nin-2,,(0)a/3 



n-1 



(0), 



2 L-z •(1;4) 

Some straightforward algebra permits to recast the result in a more compact form 

The Mellin transform defined as in d50b acts on the perturbative expression of structure function as 



(79) 



(80) 



[52 „ - s!f^]{:x.,m;z + 2n) 



where 



V, 



{0)af3 

(i;4) 



{z + 2) 



2 Do 
D 



I Vf°),f (z + 2) d^d, S^'^ (x) + Oif, Lf) 



dw 1 f d'^p l-e*"'P'^n°^(p) 



w w 



z+2 



p>m 



{2TTy p2 



pa 



(81) 



(82) 



The explicit evaluation of this integral is performed in appendix ICl 
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The relations 



x^d"^ {x^'^Yji) = [2n{2n + d-2) + d-2)]x^''Yji 

■dp{x^^Y,i) = 2nx^^Y,i 



(83) 



together with d53l yield the first order result for the leading scaling exponents for each angular component of the 
structure functions: 



C2n,j = 2n 



n{d + 2n) {d + 1) j {d + j - 2) 



d + 2 



2{d + 2){d-l) 



^+o{e) 



(84) 



The anomalous part of the scaling exponent is identified by expanding the ^ dependence of the free structure function 
upon the eddy diffusivity x^: 



P2n,j 



2n(n-l) (d + 1) j (d + j - 2) 



d + 2 



2{d + 2){d-l) 



(85) 



The result is in agreement with those given in fTOl lSHTllSll. Note that the isotropic exponent j = dominates inertial 
range scaUng. 



8.3 Second order approximation 

To second order in ^ and as Li? ^ oo we obtain the representation 



2 L 



(X) 



+e [vi^,,^''d^d, d, + i vg^gf ^^a.a, 5,9.} 5(°2(x) + o(e^ l] 

The expansion coefficients in (l86b require the evaluation of two new diagrams 



v; 



al3 

(2;4) 



• — ^wv+tAA/H— 



(86) 



v: 



af3;fi 

(2;6) 



2L>§ /■ d'^pd'^q 1 - cos[(p + q) • x] q^. n^'^(p) n°^(q) 



Z)2 y {2tt)'^^ 

q>m 
p>m 



(q + p)^ 



1. 



(87) 



Dl f d'^qd'^p (e*'!'' - l)(e"*(q+P) '' - l)(e*P '' - 1) q^ m''^(p) n°^(q) 



L»2 7 (2 7r)2'>! 

q>m 



(g2 + q . p +p2^ g2 



(88) 



Note that in tensor indices are symmetrically contracted. Consequently only the index-symmetric part of the 
above diagrams contributes to the structure function whence the omission of the semicolon between indices in (I86t . 
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The third new coefficient appearing in (l86l for L p tending to infinity is 



v; 



{2;8) 



I. 



a 11- 
(3 V. 



(89) 



Thus, the absence of momentum flow across inner scalar correlation lines reduces the evaluation of V^^g g" ^'^'^ ^'^'^ 

^(14"^ to the first order integral computed in appendix IH 

As expected by dimensional considerations, the integrals are convergent for large momentum values and logarith- 
mic for small momenta. 

The Melhn transform (l50l of (l86b acts on the individual Feynman diagrams as 



S2n - 5^°i](x,m;2; + 2n) 



^{evfJif (- + 2)+^^ [V«f (z + 2)+vg),f (z + 2)]} 



where 



,--,(0)0/3 



V, 



(2;4) 



(z + 2) 



w L>2^^+2 J (2vr)2d (q + p)2 g2 

(ir>m 



(90) 



(91) 



q>m 
p>m 



Ad'^p (e*""! '' - i)(e-™(q+p) '' _ l)(e™P '' - 1) q^ n^'^(p) n"'3(q) 



2d 



(^2 + q . p +p2) g,2 



pd 



(92) 



Guidelines for the evaluation and explicit expressions of the Mellin transform of the diagrams are deferred to ap- 
pendix [D] 

Representing ( I86t in terms of hyperspherical harmonics 15^ yields the the scaling exponents of the structure 
function for each value of the angular momentum j: 



C,2n,j =2n 



n{d + 2n) {d + I) j {d + j - 2) 
d + 2 2{d + 2){d-l) 



4(16 -5d)n'^ 



(d-l)(d + 2)3(d + 4) 



+ 



4 [4 - (j -2)j + d^ (5 + 2 j) + dj (2 j - 5) - 9d] n - 48 (d - 1) - 9 d j (d + j - 2) 



(d - 1)2 (2 + d)3 (4 + d) 



+ (n-l)Hyp2i (1,1,2 + ^,1 



3 d (6 + d) n 



+ 6n 



n(4 + d(7 + d)) -4 



(2 + d)3(4 + d)(l + d) (d2 _ 1) (2 + d)3 (4 + rf) 



3 (d3 + 6d2 + d - 4) j (d + j - 2) 



2(d- l)2(2 + d)3(4 + d) (d + 1) 
with Hyp2i denotes the Gauss hypergeometric series Q: 



HyP2i(a, ft, c, x) = > , — vFT — — p— 

r(a)r(6) ^ r(c + n)r(n + l) 



(93) 



(94) 
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The result is in agreement with those of ref.'s 10 |2l |3 derived using the ultra-violet renormalization group. From 
the computational point of view the Mellin transform applied here to the perturbative expansion does not provide the 
simplest scheme to derive the scaling exponents. It is conceptually important because it shows how to relate zero 
modes of the Hopf equations to the more general diagrammatic expansion. 

For completeness sake, it is worth noticing that gluing together truncated scalar correlation lines of dSTl . dSSl 
and ( l89t permits for any given structure function to reconstruct the diagrammatic expression of the non-universal 
contributions vanishing in the limit Lp tending to infinity. For example, the full diagrammatic expression of the 
fourth order structure function is 



54 (x; m) = 3Xx 



+3^ Txl^_Q 



• — wJ 
• — wf 



+ 6 k Xxi^^o + e' J, 



+ 3 ^ Tx|^_Q 




+12 ( e 1x1^=0 + r2'x— 



=0 / • — vvv 



+6e Xxi^=o 



+12 e Tx 



15=0 



+24 



+ 24^ 1^1=0 



+ 24^ XxU=o 




+o{e 



(95) 



9 Gradient correlations and the role of the dissipative scale 

From the short distance behavior of the structure functions of the scalar for zero molecular diffusivity we infer that 
the the field 0(x) will not be differentiable: it is only Holder continuous with exponent 1 — 0{^). Hence correlation 
functions of the gradients of 6 should blow up as the dissipative scale is taken to zero. 

The scaling Ansatz d47l imposes the existence of a precise relation between the rate of these divergences and the 
scaling exponents of the structure functions. It is convenient to illustrate the argument in the fully isotropic case, the 
generalization to anisotropy being straightforward. Using the hypothesis of inertial range universality, L and £ will 
respectively the integral and dissipative scale of the scalar field disregarding of the mechanism responsible for their 
onset. 

Radial gradient correlations at equal points and structure functions are related by the incremental ratio 



X 



a«0(x) 



2n 



y= lim -< 



^(x) - 6(0) 



2n 



lima;-2"S2n(x,^, L) 

X 10 



(96) 



By scaling 



S2n{x,e,L)=x^''r^^cn{L/e). 



Supposing 

and matching atx = £ with the scaling Ansatz ( l?7l we infer 

/ T \ P2 7 



(97) 



(98) 



17 



As the dissipative scale tends to zero gradient correlations are seen to blow up at equal points as a power law with 
exponent determined by the anomalous scaling exponent of the structure functions. Thus ( l98l relates the existence of 
anomalous scaling and the dissipative anomaly in the energy flux of the scalar field. Taking the angular average of 
( 198 1 relates it to the dissipative anomaly 



(99) 



where the angular average is given by 



4"r(n + i)r (f + 



9.1 Perturbative expansion for radial gradient correlations 



(100) 



Eq. (|98j) suggests that it should be possible to determine the scaling exponents from a perturbative expansion of 
radial gradient correlations rather than of the structure functions fSl. By dimensional analysis a similar expansion is 
seen to generate Feynman diagrams logarithmic at all momentum scales. The evaluation of the associated integrals is 
therefore greatly simplified. 
The identity 



permits to derive the perturbative expansion of radial gradients from 



(101) 



u. 



(1)q/3 



(1;4) ^""^/^ +^(2;4) '^"'^/3 



(0)a/3 



+e ^U^^^^^'-d^dpd, + \u^^.;^^U^^.;^'^^^p^,^^^ < i^-d^ef- yc +0{e) (102) 

where -< • >-g denotes averaging with respect to the Gaussian measure. The argument of the logarithms appearing 
in the expansion is now the ratio between the integral and dissipative scales of the velocity field. As in section |8]the 
result is written by keeping fixed the eddy diffusivity (fTTl and in the limit of integral scale of the forcing L p tending 
to infinity. 

The ^//-coefficients in (I102t are related through the general definition 



(n)._ 



u: 



(103) 



to the evaluation of the Feynman diagrams 



al3 

(i;4) 




/3 



D 



d'^p (p-x)2n°^(p) 



(d-i)(d + 2)e 



M>p>m 

m J 



(2vr) 



T"^(x,2) 



P' 



(104) 
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and 



(2;4) 




n2 

L»2 



M>q>m 
M>p>m 



d'^pd'^q [(p + q) • x]2 q^q^ Hf"" {p) 1^ {^1) 
(2 7r)2'i (q + p)2 g2 



(105) 



(2;6) 




D2 
L»2 



d'^'g (q • x) [(q + p) • x] (p • x) q^ n^'^(p) n"^(q) 



M>g>m 
M>p>m 



(2 7r)2'i (q2+q.p + p2)^2 



pd+S. qd+S, 



(106) 



The evaluation of the last coefficient appearing in ( 11021 does not require the evaluation of further integrals 



a f3;^u 

{2;8) 



a 



f3- 




"(1;4)"(1;4) 



(107) 



As in the case of structure function only the index-symmetric part of the diagrams above matters in the evaluation of 
(11021) whence the omission of the semicolon between indices. 

The Mellin transform proves again useful in extracting logarithms from ( 11021 . The transform is well defined since 
the condition M > m provides an infra-red cut-off: 



Q2ni'm,M;z) :-- 



m 



OC 



w 



MV 
m J 



(108) 



Evaluated at m equal one or equivalently by keeping constant in the ^-expansion, dlOSt permits to apply the 
formula (I53t to the computation of the scaling exponent. The explicit expressions of the Mellin transform of the 
Z//-coefficients together with some details about their computation are given appendix |F| Straightforward algebra then 
shows that from (I102t the 0(^^) result ( 1931 for the anomalous scaling exponent is recovered. Finally it is readily 
verified that (l93l can be related to the perturbative expansion of the structure function in the limit 



^2n(m,M) 



„2n 



lim 



™2n 



(109) 



Namely, the Z/^-coefficients in (|55Jl coincide with leading order of the Taylor expansion in the spatial increments of the 
V-coefficients in (l86l. 



10 Renormalization group methods and the validation of scaling Ansatze 

The zero mode analysis leading to the scaling Ansatz (P7l justifies the exponentiation of the logarithms encountered 
in the perturbative expansion in powers of ^. However, the concept of a zero mode is based on the existence of a 
closed hierarchy of Hopf equations for equal time correlations. This is a very special feature of the Kraichnan model 
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which is not present in more reaUstic models of fluid turbulence. It is therefore useful to look for a validation of 
scaling directly in framework of the MSR formalism. 

Such a validation is well known in the statistical field theory of dynamical and critical phenomena namely the 
Wilson formulation of the renormalization group (see for example fSF 57 1) and by its subsequent refinements (com- 
prehensive presentations can be found for example in the monographes 1 14, J9,,^ J3,^, 59 1). 

In the theory of critical phenomena scaling of correlation functions occurs at large spatial scales. They exhibit 
at the critical point a power law fall-off which is insensitive to the short distance details of the system. The univer- 
sality of the long distance or infra-red behavior with respect to the ultra-violet is formalized by direct or ultra-violet 
renormalization. In the coming subsection the main idea will be recalled with the aim of applying it to the Kraichnan 
model where universality is expected to emerge from the renormalization of infra-red rather than ultra-violet degrees 
of freedom once the dissipative scale has been set to zero. 

10.1 Wilson direct (ultra-violet) renormalization group 

We start by briefly recalling the Wilson renormalization group as applied to the study of correlation functions at large 
spatial scales. A field theory or a spin system is specified by an action functional A (or Hamiltonian) depending on 
fields, spins and the hke here collectively denoted by cf) where is a random field i;^>(x) with probability distribution 
formally given by 

V[cp\ (X e--^('^)p[(/>] (110) 

Local observables 0{x) are functions of cj) and its derivatives at the point x. In analogy to quantum field theory, they 
are often called operators. 

At the critical point correlation functions of local operators exhibit scaling 

~< 0(x)C'(y) |x - yj-^^o (HI) 

with rjo the scaling dimension of the operator O. Scaling becomes exact in the scaling limit i.e. for the random fields 

(x) = limX-'^^O (A^^x) (112) 

xio ^ ' 

whose two point function is 

^ a (x) a (y) ^« |x - yj-^'^o (113) 

Thus, the scaling operator (x) describes the long distance behavior of O (x). 

In statistical mechanics the fields (\) and hence O have an UV cutoff (e.g. the lattice spacing). has no such 
cutoff. Wilson's idea was to combine the scaling limit with a coarse graining operation so that the fields retain a fixed 
UV cutoff and the operation, called Renormalization Group, can be viewed as a on mapping probability distributions 
(or actions). This leads to a theory of the scaling dimensions r/o . 

In the simplest setup the fields (\) have a cutoff in momentum space e.g. defined by having Fourier transform 
with support for radial values of the momentum in [0, M]. Let consist of the Fourier components of ^ in the range 
[0, AM] and ^ the ones on [AM, M]: 

0(x) = (^(x) + 50(x) (114) 

Define a rescaling operation 

0(;,)(x) :=A-''°0(A-ix) (115) 
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and let 

<^'(X) = (X) (116) 

SO (\)' has momenta on [0, M]. Then we have the decomposition 

</>(x) = 0'(,/^)(x) + 5<^(x) (117) 

The expectation of O can be rewritten as 



Integrating over 50, this becomes 



where we defined the coarse-grained or renormahzed operator 

/ P[.50 ]0(0^(,/,)+.50)e-^('^a A)+^'^) 
/X)[(5<^]e"^^'^'(VA)+'5<^) 

and the renormahzed action functional 

{nxA){(p') = -In y"p[5(/.]e"-^^^'(VA)+^<^). (121) 

These functionals depend on the field (f)' which has momentum support on the original range [0, M]. TZx is called 
the renormahzation group transformation and clearly jO,\ is the derivative DTl\ of Tlx at A i.e. the hnearized RG 
transformation. 

The renormahzation group transformation TZx satisfies the relation 

7^Al7^A2 = 7^AlA2 (122) 

and defines a semi-group acting on the space of field functionals. 

The physical interpretation of RG is that by averaging over the ultra-violet degrees of freedom stored in Sc/) and 
rescaling the resulting theory coincides in the long distances with the original theory and differs only in irrelevant 
short distance properties. The limit A tending to zero should then describe universal long distance properties common 
to many action functionals A. It plays the role of the scahng limit in the Wilson formulation. The simplest limit is a 
fixed point 

A = 7^AA (123) 

and given an action A there is at most one choice of rj^ of the scaling dimension of the basic fields such that the 
renormahzation group flow of the action A converges to a fixed point. 

Scaling fields in the Wilson theory are local eigenoperators of the hnearized RG C\ at the fixed point Ai,: 

£^0(x) = A''o0(Ax). (124) 

These may often be found as follows. Suppose the hmits 

lim7^A^ = A (125) 
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and 



lim(/:AO)(A) = a. (126) 

AJ.0 

exist. Then satisfies eq. (1124b . 

We formulated the UV renormalization as the search of universal long distance properties of a theory with fixed 
UV cutoff. In field theory one is also interested in the continuum limit, i.e. the removal of the UV cutoff, M — > oo. 
Thus one considers a one parameter family of actions A^^ with UV cutoff M and possibly depending parametrically 
on M (through bare mass, coupling constant, wave function renormalization etc). One then fixes some momentum 
scale fh and splits the field as 

0(x) = 99(x) + 50(x) (127) 

where has UV cutoff m and the fluctuation part 50 momenta between fh and M. Call the result after integrating 
over 50 by e~-^™ . The limit 

lim A^ = Arh (128) 

M— >oo 

is the effective action of scale fh. This problem of continuum limit is obviously related to the previous one by trivial 
rescalings. For the limits to exist, A^ (after rescaling to say unit cutoff) has to approach a fixed point Ai, of the 
Wilson RG as M oo and then Arh (after rescaling again) will lie on the unstable manifold of Ai,. 

The Wilson idea can be applied with small changes to the time dependent correlation functions of solutions of 
stochastic (partial) differential equations. In that case the fields (/)(x, t) depend also on time and V is given by the 
MSR construction (and is not positive). The coarse graining takes place in space only whereas in time one scales. We 
write 

0(A)(x,t) := \-'^°0{\-^^,\'^H) (129) 
where the time exponent rjt has to be determined. In the simplest diffusion process rjt = —2. 



10.2 Infra-red renormalization 

In the theory of critical phenomena the concept of universality refers to the independence of the long distance prop- 
erties of correlation functions on the short distance details of the Hamiltonian which in the RG language translates to 
the fact that all the Hamiltonians in the domain of attraction of a given fixed point have the same critical exponents. 

In turbulence there is an inversion of scales. There are many mechanisms which may stir the onset of turbulence 
at large scales and universality refers to the independence of the inertial range scaling on these long distance details 
of the forcing. It has therefore been conjectured in the literature (see ll23l i32l for a general discussion and further 
references) that it should be possible to prove the irrelevance of infra-red degrees of freedom for inertial range scaling 
through the use of an inverse renormalization group. The adjective "inverse" must be understood in the sense that the 
renormalization group is constructed by averaging over fluctuations d(j) with support in the infra red. More precisely, 
Wilson's recursion scheme is implemented along the same lines of the direct one encoded in formulae (II 17t - (ll23t but 
with the following differences 

i The basic field is supported for radial values of the momentum in [m,(X)[ with m the infra-red cut-off. It 
is decomposed into a scaling field 0(i/a) ^^^^ support in [Am, oo[ and fluctuating field 6(p with support in 
[m, Xm[. In view of the inversion, the asymptotic regime is reached now for large values of A. 

ii If the (inversely) renormalized actions Ax converge to a fixed point Ai, as A tends to infinity then Ai, describes 
the universal short distance properties of the theory. 

In the following section the idea of inverse renormalization will be applied to the Kraichnan model to provide a 
validation of scaling complementary to the zero modes picture of section |5l 
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11 Inverse (infra-red) renormalization group for the Kraichnan model 



The idea to investigate the Kraichnan model by implementing a Wilson's infra-red recursion scheme was put forward 
first in ref.'s ll3ni32ll . There, it was argued that the Kraichnan model with a quasi-Lagrangian velocity field 1 12 1 had a 
inverse RG fixed point and scaling fields with dimensions given by the exponents found from the zero mode analysis 
of structure functions. 

The purpose here it to carry out such an analysis in more detail to the canonical Eulerian representation of the 
Kraichnan model. 



11.1 General settings for inverse renormalization 



It was observed in 1131] 1321 that one should consider the IRG transformation in the space of MSR actions without the 
forcing which should be treated as the correlation functions. The basic fields are 

0=(0,0",v) (130) 

The starting point is the measure 

e-^Wp [cP] = e--^^-^''''^UM0,0)dl^Di^) (131) 

where 



A,i9,9,v) :=-i{e,v''d^e)= ^- (132) 

and fi D the Gaussian measure with covariance D and fx^ is the Gaussian "measure" with "covariance" given by the 
free response function dSTt . To keep to the RG formalism discussed in Section 10 we introduce to R the same infrared 
cutoff m as in the velocity covariance D. This is only for notational convenience as can easily be checked below. 

The IRG transformation is defined by decomposing the fields as follows. The scaling exponents of time and the 
velocity are tied together by Galilean invariance, i.e. they are chosen so as to preserve the material derivative 

Vt = dt + v'^d^ (133) 

under scaling. This means that 

r]^ = -l-r]t. (134) 
The velocity field is decomposed into a scaling and fluctuating part 

V = \i/X) + ^v. (135) 
It is readily seen that v' is again Gaussian with covariance 

^ 't;'"(x, t)v"^{y, s) y= \-'^"'-^5{t - s) D"^(x, m) (136) 
whereas the covariance of the velocity fluctuation denoted by 

5 f{^,t)5vf {y,s) ■.= 5{t-s)6D"'^{x-y,m) (137) 

is given as 

SD"<'{^,m) = DoC J -^^^a^ip)xim,xm]{p)- (138) 
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Similar decomposition of the fields 9 and 9 

^ = ^(i/A) + (139a) 



G = 0[i/x) + SO (139b) 
-< 9'{x,t)e'{y,s) y= X-'i'>-''o+'^R{x - y,t - s, X^-^"" x^) = R'{x - jc,t - s). (140) 



Thus we will fix 

r]e + rig = d. (141) 

The fluctuation covariance is given by 

6R{-K-y,t-s)^5 9{-^,t)5e {y,s) =iH,{t-s) j ^ e*P(-y)-*P^ (p) (142) 

The IRG transformation of the measure (I131t is given by 

^-n,A(^)jy [^/] ^ e-'^^'(^''^"''^')(i/x,jK^', 9')diiD'{^') (143) 

where 

J\!,{e', 9'y) = - log/ e-^^'(^'^»d/i5ij(50, 69)diisD{S^r) (144) 
and 9, 9 and v have been decomposed according to ( I139al ).( n'39bt and ( I135I ). 

11.2 Infra-red renormalization group flow of the action 

Let us do first the 6v integral in (11441) - 

where 

& = ^(i/A) + SO (146) 
The second exponent on the RHS equals ai 59 = {) = 59 

(^"(i/A)5«^{i/A) ,<^^"^^"(i/A)5/3^{i/A)> = X^-^^'^'ie'dJ ,5Dlf^9'd^9'). (147) 
The large scale velocity is dominated by the constant mode of the velocity field 

JDJ^^ = (a« - l) - K^(x,m) + o (148) 

where 

pd+i n„^(p)x[o,H (^') (149) 
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Due to the unstable constant mode of the velocity field the only scaling under which the IRG stabilizes as A tends to 
infinity is the Gaussian scaling 

r]t = -2. (150) 

From (I136t we see that the covariance D' then tends to zero i.e. the velocity field disappears from the action. In that 
limit then 

e--^'.('P') = lim [e''^I<I'^'^-'fdf,sRiSe,6e). (151) 

A^oo J 

But 

Jdxeda9 = J dx9[^/^^dJ[^/^^ + j dx 568^59 (152) 

due to the disjoint supports of 9'^i/x-^ and S9 in momentum space. Thus we end up with the fixed point (up to a 
constant) 

Al = Y J ^^(^J dx9da9^ . (153) 

The action (I153t coincides with 

= - In ^ e*(e,<9.e> (j54^ 
where is a velocity field with constant covariance 

^ v^{t)v^{s) y= 6{t - s)>cj'^^ (155) 
with given by (fTTb . The only nonzero correlations of the fixed point theory are the multiple response functions 

n n ^ ^ 



^e*0,w^J>^^ (156) 

3=3=0 



and thus they coincide with the velocity averages of the same response functions in a random constant velocity field. 

Using Hopf equations ( ETl the response functions of the Kraichnan model with only two times are given as heat 
kernels 

n 

-< n 9{-Ku t)9{yi, s) y= e^"(*-^)(xi, . . . , x„, yi, . . . , y„) (157) 

i=l 

of the operators 



\ i / i<3 



and the multiple time ones are simple combinations of these with various n's. Under the scaling of space and time 
where 7]t = —2 the d"'^ will contract by L^^ thus explaining the fixed point. 

The reason for the trivial fixed point of the IRG lies, as stressed in ITTll^ . in the constant mode of the velocity 
field. As we saw in Section 4 this mode decouples in the stationary equal time correlators i.e. the second term drops 
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out from the operators (I158t when acting on translationally invariant functions. However, it doesn't decouple in the 
unequal time stationary correlators. For example, the two-point function of the scalar is given by 

-< 9{x, t)e{y, s) y= e-^i(*-") -< 0(x, s)9{y, s) >- (159) 

and the time dependence is dominated by the eddy diffusivity = D mT^. 

It is possible to modify the velocity covariance so that the zero mode decouples also from the time evolution while 
the stationary state is left unchanged 10111321 . an example being the quasi-Lagrangean velocity field. Then one may 
choose scaling r]t = — 2 + ^ and the IRG fixed point is less trivial. However, in both cases the scaling properties of the 
fixed point action bear very little information about the structure functions S2n- To understand their scaling we need 
to find the relevant scaling fields. In the next section we show how they appear in our model of Eulerian velocities 
thus bypassing the more involved quasi-Lagrangean formalism proposed in I31ll32l . 



12 Infra-red renormalization group analysis of the structure functions 

By definition structure functions are the averages 

52„(x) =^ [0(x) - em'- ^'^^uT ^ ^^^^^ 

with respect to the measure (11311) and the Gaussian measure of the forcing field /. For the renormalization group 
calculations it is convenient to take the forcing covariance ^ as 

^ /(pi,tl)/(P2,t2) ^= ^W(p^ + S(t - S) (161) 

where mp = Lp^. The forcing is in such a case isotropic. An example of anisotropic forcing is 

with \q\* = mp. 

As in sections [8] and 121 we work under the assumption 

^ « 1. (163) 

m 

It should be stressed that this is for simplicity only and the conclusions of the renormalization group analysis of 
the structure function hold true for an arbitrary value of the ratio of the integral scales provided the requirement of 
infra-red forcing is satisfied L45J . In the language of renormalization group this latter requirement means that 

f = 6f (164) 

which obviously is true as soon as mp G [0, Am]. Thus the structure functions are given as the expectation value 

S2n{-^-A) = [J{^)\ (165) 

where F is the spatial part of the forcing correlation and using invariance under x — > — x we replaced the scalar 
increment by 

J(x) = ^ (0(x, 0) + 0(-x, 0) - 2 ^(0, 0)) (166) 
In the limit mp tending to zero, the structure function operator involves of powers of 

/oo _ 
dt5e{0,t)6e{0,t) (167) 
-oo 

which is a field functional local in momentum space. 
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12.1 Scaling field for the structure functions 

At zeroth order in ^ the isotropic component of the renormahzed structure function has the scahng limit 



for 



12.1.1 First order in ^ 



(2n,j = 2 n 



(168) 



(169) 



We will work out perturbative corrections to ( I168t in the Umit mp ^ 0. In this limit there is no momentum flow 
along scalar correlation lines even at finite A . In consequence, couplings generated by Wick contractions of scalar 
fields in the interaction vertex ( 11321 ) with ghost fields in the forcing vertex ( I167t factorize as 

d'^p 2cosz(f-p)-l 



I I 

(^(i/A) 



(2vr) 



iPa- 



2p2 



-F{p) 



(170) 



which vanishes because of parity (recall that F{p) oc 6{p) as mp 0). At leading order in ^ the scaling field 
associated to the structure function operation is given by 



lim A^^"- / dl^,y;o(x) {CxS2n) (y; 

AToo J ^ 



+o{e 



(171) 



where 



d_ 

di 



{C\S2n) (x; 



if 

2 V 2n - 2 
2n 



S2n-2 (x) 



5v{°\f (x, A) + 2j|J),';^(x, <A(i/A))] {dad^S2) (x) 

+^ { 2n- 4 J Km) ^) + ^{tf '^(i/A))] {dc.S2) (x) (9^52) (x) 

The vertices are given by the small scale field functional 



5J"iJ)(x,(/)(l/A)) = J(^i/;,)(x)(6'(i/A) ,i;[\/;,))(^(l/A) ,^^fi/;,)) : = 



and the averaged value of its large scale (small momentum) counterpart 



2DQm^ f d'^p 1 - e*P-^n"^(p) 



D 



{2-kY 



p 



4H 



X]m,Xm] (P) = *Pa2x 



(172) 



(173) 



(174) 



whilst the notation dTTt is adopted for derivatives with respect to the Holder exponent ^. The symbol *Pa in < ll74t 
denotes the restriction of the momentum support of the velocity correlation. As A tends to infinity the diagram has the 
limit 

„2 



Afoo 



i'^^(d-l)(d + 2)^ (^'2)+'^K(l;4) W 



(175) 



27 



Little algebra yields the final form of (11711) 

Urn J (iOrfy/o(x) {CxS2n) (x; ^) = J dndY*o{^)S2n{^) 



(0)q/3, 



(176) 



provided C2n j is given by dS^ . 
12.1.2 Second order in ^ 

The calculation of the structure function scaling field can be inferred from the perturbation theory of section |8] The 
field independent part of the scaling field is given by the structure function diagrams in the presence of a fixed arbitrary 
ultra-violet cutoff. The field dependent part is obtained by pruning lines in these diagrams and replacing them with 
the corresponding pair of ultra-violet scaling fields. Accordingly, up to second order accuracy in ^, the renormalized 
structure function takes the form 



r=0 



+y {^^4) "(x, <Pii/x)) + (X, A)} d^d, 5fi(x) 

+e '^(x, 0(1/,)) + 5vg),f '^(x, A)} d^d, 5,5fi(x) 



+J {^K '''(^' <^{i/A)) + <°?8f^(^' ^)} 9ad, d,d. 5(°2(x) + 0{e, L-^') (177) 

The (5V*'s denote the diagrams dSTl . dSSl . and d88l where velocity correlations have momentum support in ]m , Am]. 
Infra-red logarithmic behavior of the field dependent vertices is identified by considering the scaling limits 



lim A^ 53' 

Ajoo 



.a/3 /X 



(1/A) 



•-^l(W^h^AA/H— Q: 



limA'Tx * 



Ajoo 



2(c; + l)lnA 



(178) 



and 



and 



/Li a 



lim A^^J"^;'^ 



lim A^T: 




lim A^SJ"''-^'^ 



Atoo 



^(2;8) 



A' 



(l/A) 



a /i- 



lim A Zx 

Afoo 



(179) 



(180) 



28 



The tensor structure of the diagrams is specified by 

Qf. ■■= \d,d.x^T^\^^2) = 6-%^ - ^"^'^"^If'^" (181) 

Only the index symmetric part of the above diagrams contributes to ( 11771 ) owing to the contraction with fully sym- 
metric quantities. As in sections ^ and sections ^ this fact is emphasized by omitting semicolons between non- 
symmetric indices. 

In section ^ it was shown that up to second order in ^ the perturbative expression of any structure function in the 
inertial range is compatible with that of a homogeneous function. This information together with the scaling limits 
(I178I) . ( 11791 ) dlSOt permit to verify after some straightforward algebra that the structure function operator has a finite 
scaling limit 



+^ (EtT^^) / 'if^.^o(x){5£)(x;0)- 5£)(x)C«lnA} 



2 



dndY*,{±) {^(J (x; 0) - 2 SS (x; 0)c£);, In A + S^i (x) [(c£^;, In A)^ - In a] } 
+0{e) = SiJ^;<P) (182) 



with 52*^ (x; (/>), i = 1,2 obtained by gathering the field dependence in (I177t according to its asymptotic behavior in 
A and provided the scaling dimension C,2n,j is specified by ( l93t . 

The conclusion of the above analysis is that the infra-red renormalization of the structure function operator pro- 
duces a well defined scaling field. It is worth noting that this conclusion holds also had we averaged out the velocity 
field from the outset. The resulting scaling field will then depend only on the ultra-violet degrees of freedom of the 
scalar and ghost fields but it has the same scaling exponent. 



13 Ultra-violet renormalization group 

Ultra-violet renormalization group addresses the question of removing the ultraviolet cutoff in the theory stemming 
from the one (M) in the velocity covariance. Although the correlation functions of the 9 fields have a well defined 
M — > oo limit the gradients will not have as was discussed in Section 9. Ultra-violet renormalization will study that 
divergence by finding the appropriate scaling fields. 

Ultra-violet renormalization was applied to the Kraichnan model in ref 's |5l|2j|3 in the framework of the minimal 
subtraction scheme (T9[ l40l ISTl l59l . The minimal subtraction scheme has the merit to provide probably the most 
computationally efficient setting for the determination of the scaling exponents C,2n,j which in |2l|3|| were determined 
up to third order in ^. 

The purpose of the present section is to reproduce to leading order in ^ the same calculation using Wilson's 
original scheme in order to render the comparison with infra-red renormalization more transparent. 

13.1 The effective action 

In order to inquire the limit M tending to infinity, it is more convenient not to use rescalings in the renormalization 
group, i.e. to proceed as in eq. (11141) . The fluctuation covariances have momenta on [rh, M] and they become in the 
limit of infinite ultra-violet cutoff 

hm 5i?(x,i) = ^Ho {t) J ^ e^P ''-"*^'*X[^,oo] (p) (183) 
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and 



(^^n"/'(p)X[^,oo](p) (184) 

Since the main interest is to determine the statistical properties of the scalar and the ghost fields, the renormalization 
group transformation will be applied directly to the interaction term obtained by averaging out all the degrees of 
freedom of the velocity field. The effective action at scale in stabilizes trivially as M ^ oo at all orders in perturbation 
theory. At first order it is simply given by 

Arn = + 0(C') (185) 

with 

D'^^ (186) 



13.2 Ultra-violet renormalization group analysis of radial gradients 

More interesting flow is found once we look at the scaling fields involving gradients of the scalar. Let Cm be the 
linearization of the above renormalization group in the limit as M tends to infinity. Acting on radial gradients ultra- 
violet renormalization gives 

£^^[x.5^(y,t)]2- = [x.50(y,t)]2" 

-^e ( '^l' ) [x ■ gg(y,t)]^"-^ [x • d 5e{y,t)] {69 d^e , D^^^ "''edf.d) 

+0{f) (187) 



In order to streamline the notation, on the right hand side of (11871) coarse grained fields with ultra-violet cutoff m are 
represented by the letters 6, 9. Leading order corrections are specified by the couplings 



(188) 



and 



[i • d 6e{y,t)f {sedge , w ° ^m pO) = 

2 



/ dsfl d\ Z)(o) °i (yi - y2) n • {dy6R) (y - y„ t - s)] 0(y„ s) (189) 

J —OO •_i •_i 



i=l j=l 

Let us shift the variables yj ^ yj + y and Taylor expand 



OO ^ 

9{y + yi,t + s)= [(y,.9y)"(s9i)'=^](y,t) (190) 



nlkl 

n,k=0 



30 



for i = 1,2. Then the leading order of this expansion gives the most singular part of the integral as the ultra-violet 
cut-off tends to infinity. Higher orders improve the ultra-violet behavior of the integral. Thus it is found 



[x-d69iy,t)] ((5^a«0,Z)W"^^a^ 



ds H d\ X • idy5R) (yi, s) D^' (y, - y^) + O 



i=l 



(191) 



by parity and 



[x • d5e{y,t)]^ {6edo,9,D^^^'''^ 698/39) 



Z^JJi^'^Cm, M) [8^9 (y, t)] [9^0 (y, t)] + O I ^ 



1 1 



(192) 



where was given in (11041) and it is here evaluated in the momentum range [m, M]. The terms neglected are 

irrelevant for the determination of the scaling dimension. This latter is exhibited in the scaling limit which in the 
present context can be taken by fixing the ratio 



m 

< ^ = T7 < 1 
M 



(193) 



and by considering 



lim A''^2".^ / dnYjo{i)CxM[x- 89^'' 

M'loo J 



An\\89\f^{l+^ 



■ (1) n{d + 2n) {d + l)j{d + j -2) 



d + 2 



2{d-l){d + 2) 



where 



(1) A 



C=0 



^(96)2" 



(194) 



(195) 



The derivation of ( 11941 exploits the integral identity 

dn y,, o(x) [x • 89{y, t)]'" = Anj [89 ■ 89^ (y, t) 
with An satisfying (see appendix |Elfor the proof) 



d-l 
2n- 1 



d + 2n 



id + l)j{d + j-2) 
2n{d-l) 



(196) 



(197) 



Existence of scaling requires the cancellation of the A dependence in (I194t whence it follows 

(1) _n{d + 2n) {d + l)j{d + j -2) 



d+2 



2{d-l){d + 2) 



(198) 



Thus, 
cut-off 



2n 



+ 0(^^) is a scaling field and the isotropic component of the radial gradients scales with the ultraviolet 



]\ f\ P2n,j 



m J 



(199) 
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with p2n,j given (l85b . 

Let us finish this section by comparing our treatment of the UV problem with that of ref 's |5l|2l|3l. In these papers 
the anomalous exponents of the scalar gradients are computed using the field theoretic RG derived from dimensionally 
regularized perturbation expansion. Their starting point for the MSR theory differs from our one in two ways. First 
their velocity covariance ^ has Dq ^ replaced by a constant with no explicit ^ dependence. Second, perturbation 
theory is done to the Stratonovich representation of the model ©. 

The first difference means that the "bare" correlation functions C2n need to be multiplied by a "renormalization 
constant" proportional to ^~2n ^ nontrivial limit as ^ ^ 0. 

The second difference leads to a logarithmic UV divergence at ^ = which can be traced to the M dependence 
of the effective diffusivity k in d26b . Since the authors work in dimensional regularization an UV cutoff doesn't enter 
but its role is played by ^ that plays the same role as d — 4 in dimensional regularization. Using the Ito representation 
in the perturbation as we do the logarithmically UV divergent tadpole diagram doesn't enter and the only divergences 
to be dealt with are infrared. In both approaches the UV problem for the action is trivial, indeed, the authors find a 
simple fixed point ^ for a running coupling constant describing the strength of the nonlinearity and their perturbation 
expansion becomes an expansion in powers of ^. 

14 Conclusions 

The Kraichnan model has the rare feature to allow for detailed analytical study of a turbulent system. In particular it 
allows to address the question of what kind of renormalization group if any is appropriate for turbulence, the traditional 
ultraviolet one or the more exotic inverse or infrared one. 

In the context of the Kraichnan model both direct and inverse renormalization can be successfully applied. Direct 
renormalization is natural when studying the short distance singularities that appear when the dissipative scale is taken 
to zero. The scaling fields are local operators in the derivatives of the scalar and the exponents may be computed using 
various versions of the UV RG, in dimensional regularization as in ref 's |5| |55l |2l |3l or in the Wilsonian framework 
as in the present paper. 

The inverse RG appears more natural for dealing with inertial range quantities such as the nonlocal operators that 
enter the study of the structure functions. The scaling fields are now very nonlocal in position space. We implemented 
the inverse RG in the Wilsonian framework. Although cumbersome it has the advantage of being conceptually clear. 
However other more computationally effective inverse schemes should be possible too. 

A natural question is whether infra-red renormalization may be of use in the analysis of physical models other 
than Kraichnan's. It is important here to stress that the possibility of successfully applying infra-red renormalization 
ultimately relies on the physics of the system. It remains a challenge for the future to establish whether such tool may 
prove useful to inquire scaling properties of systems for which direct renormalization cannot be applied. 
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Appendices 



A Mellin transform 

We use the following definition of the Mellin transform of a function / : [0, oo] R: 

fix,z)= —i^^=x' —i±^~x'Jiz). ©2) 
Jq WW'' Jo w w^ 

Suppose / is given by 

f{x)=x''g{x) ©3) 

with g decaying faster than any power at infinity. Then the Mellin transform is defined and analytic for Rez < a. It 
has a pole at a. Indeed, write 

, ,f°^dwg{w) X' p ^ -n / ^T^^. 

f{x;z)=x' — = / due ""gie^-^) ©4) 

Jq WW'''' z-a y_oo 

so the residue is given by 

Resfix; a) = -x^ g{0) ^5) 

If g is smooth at origin further "ultra-violet" poles will occur in a + n for any integer n, in correspondence with the 
Taylor expansion of g. If g has only power law decay at infinity, say 

c/(x)~x"^ 6>0 (|Aj6) 

the Mellin transform will exhibit "infra-red" poles also in — 6 — n, n = 0, .... 
The inverse Mellin transform is 

fix) = / -^xV(z) ©7) 

with c less than a. The inverse is usually computed by Cauchy's residue theorem. 
Finally for a function /(x) on R" we define 

Jq WW'' Jo W W" 

with X = |x|. 

More details can be found for example in the textbook fW]. Useful examples of application of the Mellin transform 
to diagrammatic expansions in field theory are given also in [481. 

B Evaluation of the inertial range asymptotics of the velocity field 

An algorithmically efficient evaluation of (fTSb is achieved using the integral representation of a power law: 

1 /•°° dnnt e""^" 



X^ Jq U 



r(f) ' 



z > ©2) 
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Namely, for negative values of z the order of integration in (fTSt can be inverted: 

The integral is absolutely convergent in the infra-red and is convergent in the ultra-violet because of the oscillatory 
exponential. Furthermore incompressibility constrains the tensorial structure of the integral to the form 

5"^(x,m;z) = ^jif:"^] T"/^(x,z) ©4) 

with T"'^ defined by (fTSb . Thus, the knowledge of the trace allows to reconstruct ( fT9l . The trace can be computed 
using the identity 



d^q e*5-^' _ r(-f) 



which is derived using (IBI^2I) . The result is 

■D\(x,m;z) _ Djm^-^x^ (z + d - 1) d T (f ) T (l - f ) 
T';.(x,z) z(z-e) + 2^r(^) 

with D now given by (fTTl . Inserting in (IB Ml) recovers ( fT^ . 

An alternative derivation of (IB 1^61 1 which does not use oscillatory integrals is available from l53l . 



(|B}5) 



(|B}6) 



C First order approximation 

Rewriting ( 1781 as the difference 



V("^f4)(x,m) 



d'^q n°^(q) ^Dom^ 



q>m 



n°'3(q) 



q>in 



(02) 



permits to evaluate it using the same method expounded in appendix |B] The Mellin transform of the second interval 
on the right hand side of iCHi is well defined for 5?z < —2. However, the residue of the pole in z = —2 cancels 
exactly with the first integral on the right hand side of (ICI^2t . Observing that 9C[2i has vanishing divergence its Mellin 
transform can be finally written as 



V"^(x,m;z + 2) 



T"^(x,z + 2), 



with 



(l;4)7(^''^;^ + 2) 



2{d + z + l)c{z) 



-2 <^z < 



m X 



2+z 



r^^(x;z + 2) 



{z + 2){d + z-l){d + 2 + z)z{z-^) 



(03) 



(04) 



and c{z) defined by d20b . Setting to zero, iC[3i has a double pole at z equal zero. The amplitude of such pole 
specifies is equal to minus the prefactor of the logarithm in the short distance asymptotics of the diagram 



{d+l)x^ T"^(x,2) 
(d-l)(d + 2) 



In 



m X 



(05) 
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with 



i/j{x) 



d_ 

du 



lnr(u + x) 



(06) 



M=0 



Differentiating the MeUin transform with respect to ^ higher orders contribution are found. In particular 



(d+ T°^(x,2) 
2{d-l){d + 2) 



In 



m X 



(d-l)(d + 2) 
For more details the reader is referred to 1531 



2(. + 2)l'"(^)+« 



(07) 



D Second order approximation 

The scope of this appendix is to expound the computational strategy followed in the evaluation of the second order 
diagrams and to give the final results used in the determination of the scaling exponents. Further details together with 
the computer packages used in the practical evaluation are available from (5T\. 

The evaluation of higher order integrals is hampered by the fact that the Mellin transform with respect of the spatial 
argument of the structure functions does not remove uniformly the mass cutoff from all the momentum integrations. 
The problem is obviated by taking the convolution of as many Mellin transforms as the number of momentum in- 
tegration affected by the cutoff. The procedure can be illustrated by considering the general form of second order 
diagrams: 

p ,q>m 

The integer n is fixed by the degree of homogeneity of the function /: 

f{wp,wq,wx) = w'' f{p,q,x) (03) 

Thus the translation of the origin in the complex z-plane makes sure that the integration contour of the MelUn anti- 
transform poles lies to the left of the pole in z equal zero. By rescaling the integral becomes 



oo 



dw f d'^pd'^q /(p, q, x) 

w'uj^ J (27r)2"' ^ 



V(x,. + n)=/ / z<0 04) 



p ,q>w 

Taking the convolution with a second Mellin transform allows to deal with unbounded momentum integrations: 

dw f dC, f°° du 1 f d'^pd'^q f{p,q,x.) 

V(x, z + n) — ' 



w J (2 in) Jq u u'^ j {2ttY'^ q^p'^ 

2<3tC<0 p,q>w 
— oo<Q(^<oo 

dC f d'^pd'^q /(p,q,x) 



(2 vrz) C{z-C) J (2vr)2'^ qd+Cpd+z-C 
zOiC<0 
-oo<3(^<oo 



(05) 



Performing, eventually with the help of the representation of a power-law (IBI^2t . the integration over momenta leaves 
with an integral over the Mellin variable (. In order to determine the scaling exponents up to second order, it is 
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necessary to compute the Mellin transform of each diagram up to 0(z ^) accuracy. The observation allows for some 
simplifications. Namely dimensional considerations impose 

d'^pd'^q /(p,q,x) 



/(x,C,^-C) : = 



(2vr)2rf qd+Cpd+z-C 



/(x,C,^-C) 



©6) 



with / not vanishing for z equal zero. Thus, it is possible to infer the form of the Laurent expansion for z and C in the 
neighborhood of zero 



/(x,C,^-C) 



c 



©7) 



The knowledge of the first two poles of / around C equal zero suffices to reconstruct the first two terms of the Laurent 
series around z equal zero of V: 



V(x, z + n) 



X [mx 



x'^imxY 



dC 



m 



-oo<3(^<oo 



2iTizC{z-C) 



c 



/(-3;l)(x) + /(-3;2)(x) /(-2)(x) 



+ 



+ Oiz- 



©8) 



In essence, inner Mellin transforms work as regularized Taylor expansions for the argument of multi-loop integrals. 
The method was applied to the evaluation of V(2;4) and V(2;6)- The integral V(2;8) in the limit of infinite integral 
scale of the forcing reduces to the product of two first order integrals. The evaluation of its Mellin transform through 
convolutions can be used to check of the procedure. 



D.l Evaluation of V(2;4) 

The diagram niost conveniently evaluated if its real space representation 

2 



^i.% = -J^ 'ds Id^y' 



(x-y)^ 

g 4 >c^ s 



m) 



©9) 



(4 TT s) 2 

is adopted as a starting point for its evaluation. The choice of the prefactor guarantees that at zero molecular viscosity 
the interaction term is order is of the order 0{^^). The reason is that in real space space (IDI^9t has the same structure 

Once the Mellin transform of this latter is known, the use of the Mellin convolution 



(i;4)- 



of the first order vertex V| 

technique (appendix [AJ reduces the evaluation (IDI^9I) to that of an integral similar of the same type of dTSb . More 
explicitly the Mellin transform is 



^(^ + 2) = 
2 f dC 

2<SRC<0 

— oo<Q^<oo 



ds / d'^y 



V 

e 4><* 



_(x-y)l 

g 4 >c~k s 



5'^^(y,m,C)(5,.a,Vj;Jp(y,m,z-C) dBlO) 



(4 vr s) 



The integrals over space-time variables can be performed exactly: 

(2;4)^ 



2-'m'~-^^^x'^{^^f^(l-(]T(^±^^ [Po{z,C,d)x'^6°^ -Pi{z,C,d)x^x^] 



2<SRC<0 
-oo<Qf <oo 



(d-irc{C-z)z{2 + z) (2 + d + z)(C-e)(C-^ + 0r(^)r 



2+d-C+z 



dBii) 
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with 

Po{z,C,d) = z[d^ {2 + z) + -3z -4-3d]+C[2 + d'^ z - z^ + d{2 + z + z^)] 



(Ol2a) 



Pi (z, C,d) = (2 + z)[d^ C - (2 + z + dCil + z)] 



(Ol2b) 



In order to determine scaling exponents with second order accuracy it is enough to evaluate (IDH U at ^ equal zero. 
To that goal the integral over ( can be performed by applying Cauchy theorem in the complex C, plane. Since the 
contour is has clockwise orientation residues must multiplied by a minus sign. Logarithmic contributions to $D[9l are 
associated to the triple and double pole in z of the Mellin transform. These latter ones are fully specified if (IDHH is 
approximated by its first residue for C equal zero. The result can be couched into the form 



vg),f(x;z + 2) 



2^2 



V, 



af3 



{2;4;1) 



with the function tf; defined by 06t . The tensor coefficients appearing in ( IDH3t are 



and 



{d + + d - 3) a;2 (^"^ - (d^ + d - 4) x°x^ 
(d- l)2(d + 2)2 

[4 - d (d^ + 4 ^2 + d - 10)] x2,5" ^ - 8 x^x^^ 



©13) 



(Ol4) 



(2;4;2)^"^- (d-l)2(d + 2)3 

The residue in z equal zero of (IDI^9I) specifies the inertial range asymptotics of the diagram at leading order in ^. 
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D.2 Evaluation of V(2;6) 

The Mellin transform of V(2;6) can be written as the sum of three terms 



with 



v5J:^(. + 3) 



v; 



{2;6;2) 



(z + 3) 



v5Jr(-+3) = E^S!^;5(-+3) 



dw DqUI^ 2 i 



d'^qd'^p 2sin(p-x) n^''(p) n°^(q) 

w L>2?i;2-2C J (g2 +q.p+p2)g2 pd+i qd+i 

dwD'^m^-'^^ f d'^qd'^p 2 sin(q • x) q^ n^''(p) n°^(q) 



^ J)2^z-2( J (2 7r)2'>' (g2 +q.p+p2)^2 pd+^ gd+i 

q>w 
p>w 

dw D^m'^''^^ f d'^qd'^p 2 sin[(q + p) • x] q^ n^'^(p) n"^(q) 



yj j;)2^z-2i J (2 7r)2rf (g2+q.p+p2)^2 pd+iqd+i 

q>w 
p>w 



©16) 



©17) 



In order to determine the scaling exponent within second order, can be set to zero in the integrands. The three Mellin 
integrals exist separately in the complex z-plane for values of z such that ^z < —2. The sum of the three integrals 
brings about the cancellations restoring the original domain of convergence of the Mellin transform of V(2;6)- 
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The explicit evaluation of the integrals is cumbersome but can be performed using some software for symbolic 
manipulations. The packages used for the evaluation are available for free download from L53J . The final result is 



~2 



1 

3 



z- 



2z2 



^2!;^ (X) + 



'^(2;6;2) 



The tensor coefficients are 



4x" x^x^^ + + [(d- l)5°^x'^ - (J^'^x^ - 5^^x"] 

(d- l)2(d + 2)2 

^{2-6-2') •= "^(1) '^"^ + ^^(2) (X^ 5°^ + x" (5^^^) + Z;(3)X" x^ x^ 



'^(2;6;1) V-^^ 



©18) 

©19) 
©20) 



where the scalar coefficients {f (j)}?^]^ are 



^(1) 

^(2) 



2(28 + 20(i- 15(i2 -8^3 -d^) _3(i2 + 5(i-4d2 -d3)Hyp2i (1,1,2 + 



d 1^ 

2' 4; 



4(d- l)2(d + 2)3(d + 4) 
2 (-4 + 4d + 3 d2 + d3) + 3 (4 + 3 d + d2) Hyp2i (l, 1, 2 + f , i) 



(3) 



4(d-l)2(d + 2)3(d + 4) 
2(6 + d + d2) + 3dHyp2i (l,l,2 + f,i) 



(d-l)2(d + 2)3 (d + 4) 
D.3 Evaluation V(2,8) 

As shown in the main text, when the integral scale of the forcing tends to infinity, this integral factorizes to 

^2,8) - "^(1,4) "^(1,4) 



©21) 



©22) 



Hence, the knowledge of the small scale asymptotics of the first order vertex V(i 4) suffices to determine the one of 
V(2,8)- Scaling exponents are conveniently evaluated using the Mellin transform of diagrams. Noting that V(i_4) has 
the form 



V = A{x) Inx + B{x) 



©23) 



with A, B some cut-off functions having finite value in zero and vanishing at infinity, the Mellin transform of (IDI^22t 
can be represented as 



V2(x, z) = X"" 



d-ue-" 



u u B2(eT) 



For z tending to zero from below one finds 



lim V2(x, z) 

zfO 



due"" 



A\e-)^ + 2B{e-)A{e-)^ + 



u ^2(0) 



A\0)^ + 2B{0)A{0)^ + 



+ 0(/) 



©24) 



©25) 



The equality entails that 

){0)a 

(2;8) 



V^^^'^'i^-z + A) 



l"mx 



1 V>(g±^)+^(i) 

z^ 2z2 



©26) 
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with 



,«(3;f.',„, 2 (<i + 1)2 !•* r"''(x, 2)T'"'(x, 2) 



(2;8;2) (d_ 1)2(^ + 2)3 

(027b) 

E Angular integrals 

Let uj denote the azimuthal angle in a given reference frame. The projection of powers of cos uj onto hyperspherical 
harmonics with zero magnetic numbers can be evaluated by considering the generating function 



n=0 

The exponential can be expanded in hyperspherical harmonics 

oo 
j=Q 

with Nj^d a normalization factor irrelevant for the present considerations. Since 

d — 2 \ 



cos,(.) ^ g !!^BesJ (j + .) y,;o(f^) ©3) 



angular integrals are just the n-th coefficient of the Taylor expansion 



whence ( 11971 ) follows. 



F Gradient expansion integrals 

The first order integral can be readily performed 



oo>p>m 



{d + l)x^ /My 



— ] r"'5(x,2) ©2) 



z(z + C)(d-l)((i + 2) V"^y 

Second order integrals can be performed by resorting to the Mellin-convolution techniques expounded in the previous 
appendix Id] Second order integrals can be always reduced to the scalar form 

^ f dC M^G(m,C,z-C) t^^, 

„(M,™;.)= J ^—i^^il^ e3) 

— oo<QC<°° 
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with 



G{m,C,z-0 



dpdq ip{p, q) 1'°^ dq 

z Ji q 



m Jm PQ P^'-'-q^ 



©4) 



and angular degrees of freedom re-absorbed into the definition of ip. 

The integrals to perform are essentially those calculated in |3| when direct renormalization and dimensional 
regularization were applied to compute the scaling exponents of the Kraichnan model. The reader interested to the 
details of the calculations is therefore referred to 01- Here the results are presented in the 



U. 



(0)af3_fM\ 



(2;4) 



+ 



+ 0(z- 



{0)af3;tM 

(2;6) 



MY 

m J 



where V(2;4;i) and V(2;6;i) were respectively defined in (IDH4t . (IDI-19t while 



and 



where the scalar coefficients {n(j) }'i=i are 



{d-lf {d + 2) 



«(1) 

«(3) 



8(l + d) + (-12-5d + 4d2 + rf3)Hyp2i(l,l,2 + f,|) 
4{d-lf{2 + df{4 + d) 
8 (1 + d)2 - 3 (4 + 3 + ^2) Hyp^^ (1, 1, 2 + f , 1) 

4(d- l)2(2 + d)3(4 + (i) 
-4-2d + 2d2 + 3rfHyp2i(l,l,2 + f,i) 
(d - 1)2 (2 + df (4 + d) 
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